Genomic and co-expression network analyses reveal candidate genes for oil accumulation based on an introgression population in Upland cotton (Gossypium hirsutum).

KEY MESSAGE: Integrated QTL mapping and WGCNA condense the potential gene regulatory network involved in oil accumulation. A glycosyl hydrolases gene (GhHSD1) for oil biosynthesis was confirmed in Arabidopsis, which will provide useful knowledge to understand the functional mechanism of oil biosynthesis in cotton. Cotton is an economical source of edible oil for the food industry. The genetic mechanism that regulates oil biosynthesis in cottonseeds is essential for the genetic enhancement of oil content (OC). To explore the functional genomics of OC, this study utilized an interspecific backcross inbred line population to dissect the quantitative trait locus (QTL) interlinked with OC. In total, nine OC QTLs were identified, four of which were novel, and each QTL explained 3.62-34.73% of the phenotypic variation of OC. The comprehensive transcript profiling of developing cottonseeds revealed 3,646 core genes differentially expressed in both inbred parents. Functional enrichment analysis determined 43 genes were annotated with oil biosynthesis processes. Implementation of weighted gene co-expression network analysis showed that 803 differential genes had a significant correlation with the OC phenotype. Further integrated analysis identified seven important genes located in OC QTLs. Of which, the GhHSD1 gene located in stable QTL qOC-Dt3-1 exhibited the highest functional linkages with the other network genes. Phylogenetic analysis showed significant evolutionary differences in the HSD1 sequences between oilseed- and starch- crops. Furthermore, the overexpression of GhHSD1 in Arabidopsis yielded almost 6.78% higher seed oil. This study not only uncovers important genetic loci for oil accumulation in cottonseed, but also provides a set of new candidate genes that potentially influence the oil biosynthesis pathway in cottonseed.

EB1C forms dimer and interacts with protein phosphatase 2A (PP2A) to regulate fiber elongation in upland cotton (Gossypium hirsutum).

Cotton is the most economically important natural fiber crop grown in more than sixty-five countries of the world. Fiber length is the main factor affecting fiber quality, but the existing main varieties are short in length and cannot suit the higher demands of the textile industry. It is necessary to discover functional genes that enable fiber length improvement in cotton through molecular breeding. In this study, overexpression of GhEB1C in Arabidopsis thaliana significantly promotes trichomes, tap roots, and root hairs elongation. The molecular regulation of GhEB1C involves its interactions with itself and GhB'ETA, and the function of GhEB1C regulation mainly depends on the two cysteine residues located at the C-terminal. In particular, the function activity of GhEB1C protein triggered with the regulation of protein phosphatase 2A, while silencing of GhEB1C in cotton significantly influenced the fiber protrusions and elongation mechanisms., Further, influenced the expression of MYB-bHLH-WD40 complex, brassinosteroids, and jasmonic acid-related genes, which showed that transcriptional regulation of GhEB1C is indispensable for cotton fiber formation and elongation processes. Our study analyzed the brief molecular mechanism of GhEB1C regulation. Further elucidated that GhEB1C can be a potential target gene to improve cotton fiber length through transgenic breeding.

Overexpression of cotton Trihelix transcription factor GhGT-3b_A04 enhances resistance to Verticillium dahliae and affects plant growth in Arabidopsis thaliana.

Verticillium wilt is a soil-borne fungal disease that severely affects cotton fiber yield and quality. Herein, a cotton Trihelix family gene, GhGT-3b_A04, was strongly induced by the fungal pathogen Verticillium dahliae. Overexpression of the gene in Arabidopsis thaliana enhanced the plant's resistance to Verticillium wilt but inhibited the growth of rosette leaves. In addition, the primary root length, root hair number, and root hair length increased in GhGT-3b_A04-overexpressing plants. The density and length of trichomes on the rosette leaves also increased. GhGT-3b_A04 localized to the nucleus, and transcriptome analysis revealed that it induced gene expression for salicylic acid synthesis and signal transduction and activated gene expression for disease resistance. The gene expression for auxin signal transduction and trichome development was reduced in GhGT-3b_A04-overexpressing plants. Our results highlight important regulatory genes for Verticillium wilt resistance and cotton fiber quality improvement. The identification of GhGT-3b_A04 and other important regulatory genes can provide crucial reference information for future research on transgenic cotton breeding.

Discovery of Novel 3,11-Bispeptide Ester Arenobufagin Derivatives with Potential in Vivo Antitumor Activity and Reduced Cardiotoxicity.

Arenobufagin, one of the bufadienolides isolated from traditional Chinese medicine Chan'su, exhibits potent antitumor activity. However, serious toxicity and small therapeutic window limits its drug development. In the present study, to our knowledge, novel 3,11-bispeptide ester arenobufagin derivatives have been firstly designed and synthesized on the base of our previous discovery of active 3-monopeptide ester derivative. The in vitro antiproliferative activity evaluation revealed that the moiety at C3 and C11 hydroxy had an important influence on cytotoxic activity and selectivity. Compound ZM350 notably inhibited tumor growth by 58.8 % at a dose 10 mg/kg in an A549 nude mice xenograft model. Therefore, compound ZM350 also presented a concentration-dependent apoptosis induction and low inhibitory effect against both hERG potassium channel and Cav1.2 calcium channel. Our study suggests that novel 3,11-bispeptide ester derivatives will be a potential benefit to further antitumor agent development of arenobufagin.

Analysis of the MIR396 gene family and the role of MIR396b in regulating fiber length in cotton.

Cotton fiber is one of the most important natural raw materials in the world textile industry. Improving fiber yield and quality has always been the main goal. MicroRNAs, as typical small noncoding RNAs, could affect fiber length during different stages of fiber development. Based on differentially expressed microRNA in the two interspecific backcross inbred lines (BILs) with a significant difference in fiber length, we identified the miR396 gene family in the two tetraploid cotton genomes and found MIR396b_D13 as the functional precursor to produce mature miR396 during the fiber elongation stage. Among 46 target genes regulated by miR396b, the GROWTH-REGULATING FACTOR 5 gene (GRF5, Gh_A10G0492) had a differential expression level in the two BILs during fiber elongation stage. The expression patterns indicated that the miR396b-GRF5 regulatory module has a critical role in fiber development. Furthermore, virus-induced gene silencing (VIGS) of miR396b significantly produced longer fiber than the wild type, and the expression level of GRF5 showed the reverse trends of the miR396b expression level. The analysis of co-expression network for the GRF5 gene suggested that a cytochrome P450 gene functions as an allene oxide synthase (Gh_D06G0089, AOS), which plays a critical role in jasmonate biosynthetic pathway. In conclusion, our results revealed that the miR396b-GRF5 module has a critical role in fiber development. These findings provide a molecular foundation for fiber quality improvement in the future.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton.

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Transcriptome analysis and identification of genes associated with oil accumulation in upland cotton.

Cotton is not only the most important fiber crop but also the fifth most important oilseed crop in the world because of its oil-rich seeds as a byproduct of fiber production. By comparative transcriptome analysis between two germplasms with diverse oil accumulation, we reveal pieces of the gene expression network involved in the process of oil synthesis in cottonseeds. Approximately, 197.16 Gb of raw data from 30 RNA sequencing samples with 3 biological replicates were generated. Comparison of the high-oil and low-oil transcriptomes enabled the identification of 7682 differentially expressed genes (DEGs). Based on gene expression profiles relevant to triacylglycerol (TAG) biosynthesis, we proposed that the Kennedy pathway (diacylglycerol acyltransferase-catalyzed diacylglycerol to TAG) is the main pathway for oil production, rather than the phospholipid diacylglycerol acyltransferase-mediated pathway. Using weighted gene co-expression network analysis, 5312 DEGs were obtained and classified into 14 co-expression modules, including the MEblack module containing 10 genes involved in lipid metabolism. Among the DEGs in the MEblack module, GhCYSD1 was identified as a potential key player in oil biosynthesis. The overexpression of GhCYSD1 in yeast resulted in increased oil content and altered fatty acid composition. This study may not only shed more light on the underlying molecular mechanism of oil accumulation in cottonseed oil, but also provide a set of new gene for potential enhancement of oil content in cottonseeds.

Quantitative Trait Locus Analysis and Identification of Candidate Genes Affecting Seed Size and Shape in an Interspecific Backcross Inbred Line Population of Gossypium hirsutum × Gossypium barbadense.

Seed size and shape are key agronomic traits affecting seedcotton yield and seed quality in cotton (Gossypium spp.). However, the genetic mechanisms that regulate the seed physical traits in cotton are largely unknown. In this study, an interspecific backcross inbred line (BIL) population of 250 BC1F7 lines, derived from the recurrent parent Upland CRI36 (Gossypium hirsutum) and Hai7124 (Gossypium barbadense), was used to investigate the genetic basis of cotton seed physical traits via quantitative trait locus (QTL) mapping and candidate gene identification. The BILs were tested in five environments, measuring eight seed size and shape-related traits, including 100-kernel weight, kernel length width and their ratio, kernel area, kernel girth, kernel diameter, and kernel roundness. Based on 7,709 single nucleotide polymorphic (SNP) markers, a total of 49 QTLs were detected and each explained 2.91-35.01% of the phenotypic variation, including nine stable QTLs mapped in at least three environments. Based on pathway enrichment, gene annotation, genome sequence, and expression analysis, five genes encoding starch synthase 4, transcription factor PIF7 and MYC4, ubiquitin-conjugating enzyme E27, and THO complex subunit 4A were identified as candidate genes that might be associated with seed size and shape. Our research provides valuable information to improve seed physical traits in cotton breeding.

Mapping of dynamic QTLs for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in a backcross inbred line population of Upland cotton.

KEY MESSAGE: A backcross inbred line population of cotton was evaluated for Fusarium wilt race 4 resistance at different days after inoculation (DAI). Both constitutively expressed and developmentally regulated QTLs were detected. The soil-borne fungus Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) causes Fusarium wilt including seedling mortality in cotton. A backcross inbred line (BIL) population of 181 lines, derived from a bi-parental cross of moderately resistant non-recurrent Hai 7124 (Gossypium barbadense) and recurrent parent CCRI 36 (G. hirsutum), was evaluated under temperature-controlled conditions for FOV4 resistance with artificial inoculations. Based on three replicated tests evaluated at 7, 14, 21, and 28 days after inoculation (DAI), only 2-5 BILs showed lower disease severity ratings (DSR) than the parents while 22-50 BILs were more susceptible, indicating transgressive segregation toward susceptibility. Although DSR were overall congruent between DAI, there were many BILs displaying different responses to FOV4 across DAI. Genetic mapping using 7709 SNP markers identified 42 unique QTLs for four evaluation parameters- disease incidence (DI), DSR, mortality rate (MR), and area under disease progress curve (AUDPC), including 26 for two or more parameters. All five QTLs for AUDPC were co-localized with QTLs for DI, DSR, and/or MR at one or two DAI, indicating the unnecessary use of AUDPC in QTL mapping for FOV4 resistance. Those common QTLs explained the significant positive associations between parameters observed. Ten common QTLs with negative or positive additive effects were detected between DAI. DAI-specific and consistent QTLs were detected between DAI in cotton for the first time, suggesting the existence of both constitutively expressed and developmentally regulated QTLs for FOV4 resistance and the importance of evaluating genetic populations for FOV4 resistance at different growth stages.

Quantitative Trait Locus Analysis and Identification of Candidate Genes for Micronaire in an Interspecific Backcross Inbred Line Population of Gossypium hirsutum × Gossypium barbadense.

Cotton is the most important fiber crop and provides indispensable natural fibers for the textile industry. Micronaire (MIC) is determined by fiber fineness and maturity and is an important component of fiber quality. Gossypium barbadense L. possesses long, strong and fine fibers, while upland cotton (Gossypium hirsutum L.) is high yielding with high MIC and widely cultivated worldwide. To identify quantitative trait loci (QTLs) and candidate genes for MIC in G. barbadense, a population of 250 backcross inbred lines (BILs), developed from an interspecific cross of upland cotton CRI36 × Egyptian cotton (G. barbadense) Hai7124, was evaluated in 9 replicated field tests. Based on a high-density genetic map with 7709 genotyping-by-sequencing (GBS)-based single-nucleotide polymorphism (SNP) markers, 25 MIC QTLs were identified, including 12 previously described QTLs and 13 new QTLs. Importantly, two stable MIC QTLs (qMIC-D03-2 on D03 and qMIC-D08-1 on D08) were identified. Of a total of 338 genes identified within the two QTL regions, eight candidate genes with differential expression between TM-1 and Hai7124 were identified. Our research provides valuable information for improving MIC in cotton breeding.

Discovery of 3-peptide substituted arenobufagin derivatives as potent antitumor agents with low cardiotoxicity.

Active natural productscan be valuable lead compounds and numerous drugs derived from natural products have successfully entered the clinic. Arenobufagin, one of the important active components of toad venom, indicates significant antitumor activities with limited preclinical development for its strong cardiotoxicity. Ten 3-monopeptide substituted arenobufagin derivatives have been designed and synthesized. Antitumor activity and cardiotoxicity assays lead to the discovery of compound ZM226 as a potent antitumor agent with low cardiotoxicity. These findings suggest optimization of arenobufagin on position 3 maybe an efficacious strategy for the development of antitumor drug candidates derived from arenobufagin.

Analysis of the MIR160 gene family and the role of MIR160a_A05 in regulating fiber length in cotton.

MAIN CONCLUSION: The MIR160 family in Gossypium hirsutum and G. barbadense was characterized, and miR160a_A05 was found to increase cotton-fiber length by downregulating its target gene (ARF17) and several GH3 genes. Cotton fiber is the most important raw material for the textile industry. MicroRNAs are involved in regulating cotton-fiber development, but a role in fiber elongation has not been demonstrated. In this study, miR160a was found to be differentially expressed in elongating fibers between two interspecific (between Gossypium hirsutum and G. barbadense) backcross inbred lines (BILs) with different fiber lengths. The gene MIR160 colocalized with a previously mapped fiber-length quantitative trait locus. Its target gene ARF17 was differentially expressed between the two BILs during fiber elongation, but in the inverse fashion. Bioinformatics was used to analyze the MIR160 family in both G. hirsutum and G. barbadense. Moreover, qRT-PCR analysis identified MIR160a as the functional MIR160 gene encoding the miR160a precursor during fiber elongation. Using virus-induced gene silencing and overexpression, overexpressed MIR160a_A05 resulted in significantly longer fibers compared with wild type, whereas suppression of miR160 resulted in significantly shorter fibers. Expression levels of the target gene auxin-response factor 17 (ARF17) and related genes GH3 in the two BILs and/or the virus-infected plants demonstrated similar changes in response to modulation of miR160a level. Finally, overexpression or suppression of miR160 increased or decreased, respectively, the cellular level of indole-3-acetic acid, which is involved in fiber elongation. These results describe a specific regulatory mechanism for fiber elongation in cotton that can be utilized for future crop improvement.

A Comparative Genome-Wide Analysis of the R2R3-MYB Gene Family Among Four Gossypium Species and Their Sequence Variation and Association With Fiber Quality Traits in an Interspecific G. hirsutum × G. barbadense Population.

Cotton (Gossypium spp.) is the most important natural fiber crop in the world. The R2R3-MYB gene family is a large gene family involved in many plant functions including cotton fiber development. Although previous studies have reported its phylogenetic relationships, gene structures, and expression patterns in tetraploid G. hirsutum and diploid G. raimondii, little is known about the sequence variation of the members between G. hirsutum and G. barbadense and their involvement in the natural quantitative variation in fiber quality and yield. In this study, a comprehensive genome-wide comparative analysis was performed among the four Gossypium species using whole genome sequences, i.e., tetraploid G. hirsutum (AD1) and G. barbadense (AD2) as well as their likely ancestral diploid extants G. raimondii (D5) and G. arboreum (A2), leading to the identification of 406, 393, 216, and 213 R2R3-MYB genes, respectively. To elucidate whether the R2R3-MYB genes are genetically associated with fiber quality traits, 86 R2R3-MYB genes were co-localized with quantitative trait loci (QTL) hotspots for fiber quality and yield, including 42 genes localized within the fiber length QTL hotspots, in interspecific G. hirsutum × G. barbadense populations. There were 20 interspecific nonsynonymous single-nucleotide polymorphism (SNP) sites between the two tetraploid cultivated species, of which 16 developed from 11 R2R3-MYB genes were significantly correlated with fiber quality and yield in a backcross inbred population (BIL) of G. hirsutum × G. barbadense in at least one of the four field tests. Taken together, these results indicate that the sequence variation in these 11 R2R3-MYB genes is associated with the natural variation (i.e., QTL) in fiber quality and yield. Moreover, the functional SNPs of five R2R3-MYB allele pairs from the AD1 and AD2 genomes were significantly correlated with the gene expression related to fiber quality in fiber development. The results will be useful in further elucidating the role of the R2R3-MYB genes during fiber development.

Genome-wide association study of the oil content in upland cotton (Gossypium hirsutum L.) and identification of GhPRXR1, a candidate gene for a stable QTLqOC-Dt5-1.

Cottonseed oil is one of the most important renewable resources for edible oil and biodiesel. To detect QTLs associated with cottonseed oil content (OC) and identify candidate genes that regulate oil biosynthesis, a panel of upland cotton germplasm lines was selected among those previously used to perform GWASs in China. In the present study, 13 QTLs associated with 53 common SNPs on 13 chromosomes were identified in multiple environments based on 15,369 polymorphic SNPs using the Cotton63 KSNP array. Of these, the OC QTL qOC-Dt5-1 delineated by nine SNPs occurred in a confidence interval of 4 SSRs with previously reported OC QTLs. A combined transcriptome and qRT-PCR analysis revealed that a peroxidase gene (GhPRXR1) was predominantly expressed during the middle-late stage (20-35 days post anthesis) of ovule development. The overexpression of GhPRXR1 in yeast significantly increased the OC by 20.01-37.25 %. Suppression of GhPRXR1 gene expression in the virus-induced gene-silenced cotton reduced the OC by 18.11%. Our results contribute to identifying more OC QTLs and verifying a candidate gene that influences cottonseed oil biosynthesis.

QTL analysis and candidate gene identification for plant height in cotton based on an interspecific backcross inbred line population of Gossypium hirsutum × Gossypium barbadense.

KEY MESSAGE: We constructed the first high-quality and high-density genetic linkage map for an interspecific BIL population in cotton by specific-locus amplified fragment sequencing for QTL mapping. A novel gene GhPIN3 for plant height was identified in cotton. Ideal plant height (PH) is important for improving lint yield and mechanized harvesting in cotton. Most published genetic studies on cotton have focused on fibre yield and quality traits rather than PH. To facilitate the understanding of the genetic basis in PH, an interspecific backcross inbred line (BIL) population of 250 lines derived from upland cotton (Gossypium hirsutum L.) CRI36 and Egyptian cotton (G. barbadense L.) Hai7124 was used to construct a high-density genetic linkage map for quantitative trait locus (QTL) mapping. The high-density genetic map harboured 7,709 genotyping-by-sequencing (GBS)-based single nucleotide polymorphism (SNP) markers that covered 3,433.24 cM with a mean marker interval of 0.67 cM. In total, ten PH QTLs were identified and each explained 4.27-14.92% of the phenotypic variation, four of which were stable as they were mapped in at least two tests or based on best linear unbiased prediction in seven field tests. Based on functional annotation of orthologues in Arabidopsis and transcriptome data for the genes within the stable QTL regions, GhPIN3 encoding for the hormone auxin efflux carrier protein was identified as a candidate gene located in the stable QTL qPH-Dt1-1 region. A qRT-PCR analysis showed that the expression level of GhPIN3 in apical tissues was significantly higher in four short-statured cotton genotypes than that in four tall-statured cotton genotypes. Virus-induced gene silencing cotton has significantly increased PH when the expression of the GhPIN3 gene was suppressed.

Genome-Wide Identification, Sequence Variation, and Expression of the Glycerol-3-Phosphate Acyltransferase (GPAT) Gene Family in Gossypium.

Cotton is an economically important crop grown for natural fiber and seed oil production. Cottonseed oil ranks third after soybean oil and colza oil in terms of edible oilseed tonnage worldwide. Glycerol-3-phosphate acyltransferase (GPAT) genes encode enzymes involved in triacylglycerol biosynthesis in plants. In the present study, 85 predicted GPAT genes were identified from the published genome data in Gossypium. Among them, 14, 16, 28, and 27 GPAT homologs were identified in G. raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. Phylogenetic analysis revealed that a total of 108 GPAT genes from cotton, Arabidopsis and cacao could be classified into three groups. Furthermore, through comparison, the gene structure analyses indicated that GPAT genes from the same group were highly conserved between Arabidopsis and cotton. Segmental duplication could be the major driver for GPAT gene family expansion in the four cotton species above. Expression patterns of GhGPAT genes were diverse in different tissues. Most GhGPAT genes were induced or suppressed after salt or cold stress in Upland cotton. Eight GhGPAT genes were co-localized with oil and protein quantitative trait locus (QTL) regions. Thirty-two single nucleotide polymorphisms (SNPs) were detected from 12 GhGPAT genes, sixteen of which in nine GhGPAT genes were classified as synonymous, and sixteen SNPs in ten GhGPAT genes non-synonymous. Two SNP markers of the GhGPAT16 and GhGPAT26 genes were significantly correlated with cotton oil content in one of the three field tests. This study shed lights on the molecular evolutionary properties of GPAT genes in cotton, and provided reference for improvement of cotton response to abiotic stress and the genetic improvement of cotton oil content.

Genetic variation of dynamic fiber elongation and developmental quantitative trait locus mapping of fiber length in upland cotton (Gossypium hirsutum L.).

BACKGROUND: In upland cotton (Gossypium hirsutum L.), genotypes with the same mature fiber length (FL) might possess different genes and exhibit differential expression of genes related to fiber elongation at different fiber developmental stages. However, there is a lack of information on the genetic variation influencing fiber length and its quantitative trait loci (QTLs) during the fiber elongation stage. In this study, a subset of upland cotton accessions was selected based on a previous GWAS conducted in China and grown in multiple environments to determine the dynamic fiber length at 10, 15, 20, and 25 days post-anthesis (DPA) and maturity. The germplasm lines were genotyped with the Cotton 63 K Illumina single-nucleotide polymorphism (SNP) array for GWAS. RESULTS: A total of 25, 38, 57, 89 and 88 SNPs showed significant correlations with fiber length at 10, 15, 20 and 25 DPA and maturity, respectively. In addition, 60 more promising SNPs were detected in at least two tests and two FL developmental time points, and 20 SNPs were located within the confidence intervals of QTLs identified in previous studies. The fastest fiber-length growth rates were obtained at 10 to 15 DPA in 69 upland cotton lines and at 15 to 20 DPA in 14 upland cotton accessions, and 10 SNPs showed significant correlations with the fiber-length growth rate. A combined transcriptome and qRT-PCR analysis revealed that two genes (D10G1008 and D13G2037) showed differential expression between two long-fiber genotypes and two short-fiber genotypes. CONCLUSIONS: This study provides important new insights into the genetic basis of the time-dependent fiber-length trait and reveals candidate SNPs and genes for improving fiber length in upland cotton.

Genome-Scale Analysis of the WRI-Like Family in Gossypium and Functional Characterization of GhWRI1a Controlling Triacylglycerol Content.

Cotton (Gossypium spp.) is the most important natural fiber crop and the source of cottonseed oil, a basic by-product after ginning. AtWRI1 and its orthologs in several other crop species have been previously used to increase triacylglycerols in seeds and vegetative tissues. In the present study, we identified 22, 17, 9, and 11 WRI-like genes in G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, respectively. This gene family was divided into four subgroups, and a more WRI2-like subfamily was identified compared with dicotyledonous Arabidopsis. An analysis of chromosomal distributions revealed that the 22 GhWRI genes were distributed on eight At and eight Dt subgenome chromosomes. Moreover, GhWRI1a was highly expressed in ovules 20-35 days after anthesis and was selected for further functional analysis. Ectopic expression of GhWRI1a rescued the seed phenotype of a wri1-7 mutant and increased the oil content of Arabidopsis seeds. Our comprehensive genome-wide analysis of the cotton WRI-like gene family lays a solid foundation for further studies.

A genome-wide analysis of the lysophosphatidate acyltransferase (LPAAT) gene family in cotton: organization, expression, sequence variation, and association with seed oil content and fiber quality.

BACKGROUND: Lysophosphatidic acid acyltransferase (LPAAT) encoded by a multigene family is a rate-limiting enzyme in the Kennedy pathway in higher plants. Cotton is the most important natural fiber crop and one of the most important oilseed crops. However, little is known on genes coding for LPAATs involved in oil biosynthesis with regard to its genome organization, diversity, expression, natural genetic variation, and association with fiber development and oil content in cotton. RESULTS: In this study, a comprehensive genome-wide analysis in four Gossypium species with genome sequences, i.e., tetraploid G. hirsutum- AD1 and G. barbadense- AD2 and its possible ancestral diploids G. raimondii- D5 and G. arboreum- A2, identified 13, 10, 8, and 9 LPAAT genes, respectively, that were divided into four subfamilies. RNA-seq analyses of the LPAAT genes in the widely grown G. hirsutum suggest their differential expression at the transcriptional level in developing cottonseeds and fibers. Although 10 LPAAT genes were co-localised with quantitative trait loci (QTL) for cottonseed oil or protein content within a 25-cM region, only one single strand conformation polymorphic (SSCP) marker developed from a synonymous single nucleotide polymorphism (SNP) of the At-Gh13LPAAT5 gene was significantly correlated with cottonseed oil and protein contents in one of the three field tests. Moreover, transformed yeasts using the At-Gh13LPAAT5 gene with the two sequences for the SNP led to similar results, i.e., a 25-31% increase in palmitic acid and oleic acid, and a 16-29% increase in total triacylglycerol (TAG). CONCLUSIONS: The results in this study demonstrated that the natural variation in the LPAAT genes to improving cottonseed oil content and fiber quality is limited; therefore, traditional cross breeding should not expect much progress in improving cottonseed oil content or fiber quality through a marker-assisted selection for the LPAAT genes. However, enhancing the expression of one of the LPAAT genes such as At-Gh13LPAAT5 can significantly increase the production of total TAG and other fatty acids, providing an incentive for further studies into the use of LPAAT genes to increase cottonseed oil content through biotechnology.

[Genome-wide analysis of the LPAAT gene family in Gossypium raimondii and G. arboreum, and expression analysis of its orthologs in G. hirsutum].

Lysophosphatidic acid acyltransferase (LPAAT) which converts lysophosphatidic acid into phosphatidic acid is a key enzyme in biosynthesis pathway of lipid in plants. In this study, we identified 17 members of the LPAAT gene family from genomic data of G. raimondii-D5 and G. arboreum-A2. Analysis of gene structure, chromosome distribution and phylogenetic evolution of LPAAT genes in diploid Gossypium using bioinformatics approaches showed that these genes can be divided into distinct subfamilies based on the distance of their genetic relationship. Moreover, the gene structures were similar within LPAAT subfamily members. The amino acid sequences encoded by LPAAT family genes contained three conserved motifs, including ΦFPEGTR-G binding site and Φ-NHQS- ΦDΦΦ catalytic site. Phylogenetic analysis of LPAAT gene family demonstrated significant differences in evolution of LPAAT in different species. Finally, expression analysis of G. hirsutum ovules in different stages from RNA-seq and qRT-PCR data indicated that LPAAT gene may play a positive role in oil accumulation. Our studies facilitate understanding of the function of LPAAT gene family in Gossypium and selecting better LPAAT genes for further functional validation.